NSF-ITP-95-163 

A Diagrammatic Theory of Random Scattering Matrices 
for Normal— Superconducting Mesoscopic Junctions 

N. Argaman and A. Zee 
Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106, USA 

(Submitted: December 28, 1995) 

Abstract 

The planar-diagrammatic technique of large-iV random matrices is ex- 
tended to evaluate averages over the circular ensemble of unitary matrices. 
It is then applied to study transport through a disordered metallic "grain", 
attached through ideal leads to a normal electrode and to a superconducting 
electrode. The latter enforces boundary conditions which coherently couple 
electrons and holes at the Fermi energy through Andreev scattering. Con- 
sequently, the leading order of the conductance is altered, and thus changes 
much larger than e 2 /h are observed when, e.g., a weak magnetic field is ap- 
plied. This is in agreement with existing theories. The approach developed 
here is intermediate between the theory of dirty superconductors (the Usadel 
equations) and the random-matrix approach involving transmission eigenval- 
ues (e.g. the DMPK equation) in the following sense: even though one starts 
from a scattering formalism, a quantity analogous to the superconducting 
order-parameter within the system naturally arises. The method can be ap- 
plied to a variety of mesoscopic normal-superconducting structures, but for 

brevity we consider here only the case of a simple disordered N-S junction. 
PACS numbers: 73.23Ps, 74.50.+r 
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I. INTRODUCTION 



Dissipative single-electron transport in mesoscopic structures involving normal and su- 
perconducting elements has in recent years become a topic of lively experimental |l]-|6[] and 
theoretical 0-l2| activity. From both of these points of view, this is a natural extension of 
the interest in usual (i.e., normal) mesoscopic systems. Indeed, many of the technological 
and mathematical techniques used here are familiar from that field: lithography, lock-in am- 
plifiers, diagrammatic perturbation theory, and random matrices, to name a few. However, 
this extension has often lead to surprises. 

Consider for example the effect of a magnetic field on a normal mesoscopic metallic 
system, which can be described in the leading-order approximation by classical dynamics 
of non-interacting electrons at the Fermi surface. The conductance is then determined by 
the density of states and the geometric details of the sample, and is of the order of N e 2 /h 
where N is the number of propagating modes in the system, e is the electronic charge and h 
is Planck's constant. In the present work, N is taken to be large, whereas additional factors 
which appear in the conductance, such as the factor l/L where / is the transport mean-free- 
path and L is the length of the system, are taken to be of order unity and considered as 
fixed geometric characteristics of the system. 

If the magnetic field in such a system is so weak that it does not affect the classical 
motion of the electrons, it may still shift the phases of different interfering partial waves in 
the Schrodinger equation, thus leading to changes in the conductance which are well-known 
P~3|l to be universal, of the order of the quantum unit of conductance, e 2 /h. This contrasts 



sharply with the situation in the presence of one or more superconducting electrodes (we use 
"electrode" rather than the more technical nomenclature — "particle reservoir"). By the 
process of Andreev reflection Jl4|,[15|, an electron-like excitation in the normal part of the 
system may reflect off the normal-superconducting (N-S) boundary and give rise to a hole- 
like excitation. The extra charge of two electrons is thus transformed into a Cooper pair, 
which is carried away by the superconductor. The Andreev-reflected hole is necessarily found 
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to be in the state time-reversed to the impinging electron. If perfect coherence is maintained 
between electrons and holes, an analogy may be drawn between the N-S boundary and a 
phase-conjugating mirror, and one finds that the hole retraces the possibly complicated 
motion of the electron. A weak magnetic field may break the symmetry between electrons 
and holes, leading to large effects in the conductance, of the order of N e 2 /h. 

Our purpose here is to suggest a new theoretical technique for calculating such effects 
in dirty mesoscopic N-S structures, which complements the existing theoretical tools in 
interesting ways. In this approach the system is described as a collection of ergodic scatterers, 
each possesing a random scattering matrix. The scattering cavities are connected to each 
other and to the external electrodes through leads of various widths. In this first paper, 
we demonstrate the new approach by studying a system described by only one scattering 
cavity, attached to a single normal electrode and a single superconducting electrode through 
two ideal leads (see Fig. 1). The discussion of more complicated structures (having more 
than one scattering cavity, more than two electrodes, and non-ideal leads) is deferred to a 
future publication ]IB . 



At first sight, it might seem that the conductances of such systems would always be larger 
in the case with electron-hole symmetry than in the non-symmetric case, because whenever 
an electron happens to undergo Andreev reflection, the charge transport is "automatically" 
doubled by the corresponding motion of the hole. However, only some of the partial waves 
representing an electron entering the system from a normal electrode will end up impinging 
on a superconducting electrode, and thus the hole only imperfectly reproduces the wavefronts 
and the motion of the electron. Furthermore, the presence of Andreev scattering actually 
reduces the density of states at the Fermi surface — an aspect of the proximity effect which 
allows the gap of the superconductor to be manifest in the adjacent normal material — and 
thus the conductance may sometimes be lowered more than it is enhanced. In the approach 
to be described below, the magnitude of both of these effects is described by a single quantity 
/, which in principle may depend on the position within the (normal part of the) system, 
and on the energy of the pertinent electrons. We define / as the (averaged, leading-order 



of the) probability amplitude for an electron at that position and energy to propagate to 
a superconducting electrode, be Andreev-reflected, and then to coherently propagate back 
to the same position as a hole. The density of states at the Fermi level is then equal to 
(1 — |/| 2 )/(1 + |/| 2 ) times the normal one jD]] (i.e. the density of states in the absence of 
induced superconductivity). 

The present theoretical approach creates a link between the available calculational tools 
for disordered mesoscopic N-S structures. On the one hand, the starting point is the expres- 
sion of the conductance in terms of scattering properties (an appropriate Landauer formula), 
which themselves are assumed to be described well by certain ensembles of scattering ma- 
trices. In this respect, the approach is similar in spirit to a large number of recent articles 
(e.g. those studying the Dorokhov-Mello-Pereyra-Kumar equation P-|ll,18|), in which the 



distribution of transmission-matrix eigenvalues for various structures was studied and used 
to obtain the conductance |1| (the transmission matrix is simply related to the scattering 
matrix of a structure, and carries the same information). On the other hand, the large- 
N diagrams used naturally lead to the definition of the probability-amplitude /, which is 
closely related to the space- and frequency-dependent superconducting order parameter of 



the Usadel equations |2(J . This set of equations, which is widely used in the theory of dirty 
superconductors, is ultimately based on the Keldysh technique of diagrammatic perturbation 
theory. 

Another approach which proves useful for disordered mesoscopic systems is the semi- 



classical approach of the Gutzwiller formula or the van Vleck formula However, the 
semiclassical approximation for disordered systems is not strictly unitary, and this turns 
out to be a fatal flaw in the present context (see the next section). Additional approaches 
developed for mesoscopic systems in general, such as the nonlinear cx-model ||22|| , do not 
suffer from this flaw and are, in principle, applicable. 

The outline of the paper is as follows. The more technical part of the introduction is given 
in the following subsection. In section II, we develop the formalism of planar diagrams for 
averaging over large- N random unitary matrices: the diagrammatic elements are discussed, 



and used to form "self-energies" and "effective couplings". This section is general, in the 
sense that its results are used for both the simple junction studied in section III and the 



more complicated structures (networks of cavities) to be discussed in Ref . (16| . In section III 
we proceed to apply the formalism to the simplest possible system — that having only two 
leads — and compare our results with the literature. Finally, conclusions and an outlook 
are given in section IV. 



Preliminaries 

As discussed above, only normal-superconducting systems with simple geometries - 
those describable by a single "grain" (Fig. 1) — are considered in the present work. The 
grain is taken to be an ideal scattering cavity with ergodic single-electron internal dynamics, 
i.e., its scattering matrix S is assumed to be a random member of the Circular Orthogonal 
Ensemble (COE) of random matrix theory [f23] (the ensemble of unitary symmetric matrices 
- the symmetry or "orthogonality" corresponds physically to time-reversal symmetry). 
This ensemble is known to describe a wide universality class of physical systems, including 
at least two different relevant categories: clean cavities with shapes which lead to chaotic 
classical dynamics, and disordered materials for which the chaotic scattering is provided 
by impurities. In both of these cases one assumes that the escape time for an electron in 
the cavity is much longer than the time required for ergodicity — otherwise the description 
using a single random scattering matrix becomes unjustifiable. 

The grain is connected to external electrodes via leads which are assumed to be ideal 
("waveguides" of electrons), and can have different widths, denoted W n . The physical widths 
are assumed to be much larger than the Fermi wavelength, and the W n 's are thus large 
integers which give the number of propagating modes in each lead. The total number of 
modes, N = J2n W n , gives the size of the scattering matrix S. In the simplest case, to be 
discussed in detail in section III, there are only two leads — one "superconducting lead" 
of width W\ connects the system to a superconducting electrode, and one "normal lead" 
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of width W2 connects it to a normal electrode. In this case, the S matrix of the cavity is 
written as 

(r t'\ 

where the reflection matrices r (of size W\ x W\) and r' (of size W2 x W2) are not equal in size, 
and the two transmission matrices t and t' are rectangular (except for the case W\ = W2). 
In general, S may have more than four blocks, describing scattering from and into more than 
two leads (at least one normal lead and at least one superconducting lead are necessary for 
the conductances calculated here to be of relevance). However, the statistical properties of 
the ^-matrix elements, according to the COE, do not depend on the number and sizes of 
the blocks. 

Once an electron has scattered in the cavity it may enter the superconducting lead 
and impinge upon the superconducting electrode. As we will assume that the electron is 
propagating at the Fermi level, it can not enter the superconductor because of the gap, and 
it will be completely reflected. However, it may in principle undergo normal reflection, and 
return as an electron, or Andreev reflection [TJ|] and return as a hole. In the present work 



we will assume for simplicity that the probability of normal reflection vanishes, so that the 
probability amplitude for Andreev reflection is of magnitude unity (this assumption may 
be relaxed by introducing tunnel barriers in the ideal leads). The phase of this probability 
amplitude is determined from the Bogoliubov - de Gennes equations p4| , and the amplitude 



is found to be ie~ l< ^ where <fi is the phase of the superconducting order parameter in the 
electrode. The scattering of holes in the cavity is described by S*, the complex conjugate of 
the scattering matrix for electrons. The amplitude for a hole impinging on the N-S boundary 
to be Andreev-reflected back into an electron is ie 1 ^. 

The extra phase factors of i are crucial to the description of the induced proximity effects 
- they lead to a relative minus sign, i.e. destructive interference, between the amplitude for 
an electron to reach any final state through a certain history, and the amplitude to reach 
the same final state through a history with two additional Andreev reflections (into a hole 



and back into an electron). In the absence of any time-reversal symmetry-breaking effects, 
the dynamics is just such that two consecutive Andreev reflections tend to reproduce the 
same state, and so this destructive interference leads to a suppression of (again) the leading 
order expression for the density of states, by the factor of (1 — |/| 2 )/(1 + |/| 2 ) mentioned 
above. 

According to the appropriate Landauer formula ||25|| , the average conductance of such a 



structure (at zero temperature and zero bias voltage) is given by 

(vns) = ^ (trT eh T} h ) , (2) 

where T e h is an electron-hole "transmission" matrix, giving the probability amplitudes for 
an electron injected into the cavity through any one of the modes of the normal lead to come 
back as a hole in any other of the modes of that lead (the notation (. . .) denotes averaging 
of the S matrix over the ensemble). This is similar to the standard Landauer formula, 
o~n = (2e 2 /h)trt t* where the factor of e 2 /h is the quantum unit of conductance, i.e., the 
conductance per mode of an ideal lead with no scattering, the factor of 2 comes from spin 
(which will play no role in this paper), and the trace sums the probabilities for electrons to 
be transmitted through the system. The extra factor of 2 in Eq. (^j) is due to the fact that 
for every electron from the normal lead that comes back as a hole, by charge conservation 
two electrons have gone into the superconductor. 

The amplitudes in T e h are given by a sum over all possible histories. In the simplest kind 
of history, an electron injected into the cavity traverses the cavity with the transmission 
amplitude t', Andreev reflects off the superconductor with the amplitude % and becomes a 
hole, and then the hole traverses the cavity with transmission amplitude t*. This gives a 
contribution of t*it' . As the system under discussion has only a single superconducting lead, 
we have taken here for simplicity = 0, and avoided the extra phase factors (because of gauge 
invariance, they are only important in a system with at least two superconducting leads). 
In the next simplest history, after Andreev reflection, the hole is scattered back towards the 
superconductor (rather than traversing the cavity), with the reflection amplitude r*. It then 



Andreev reflects and becomes an electron. The electron bounces off the cavity back towards 
the superconductor again with the reflection amplitude r, only to come back as a hole after 
Andreev reflection. The hole finally traverses the cavity with transmission amplitude t*. 
This contributes an amplitude of t*irir*it' to T e ^. The matrix T e h may thus be written as 

T e h — t*(i + irir*i + irir*irir*i + . . .)t' (3) 

This series for T e h could be written directly in terms of the S matrix — for example, the 
term t*irir*it' can be written as i 3 S*SS*S — provided that we adopt the peculiar convention 
that in this product of four unitary matrices the internal indices are to be summed over only 
the W\ indices pertaining to the lead connecting the cavity to the superconductor and the 
external indices are restricted to the W2 indices pertaining to the external lead. 

Several simplifying assumptions are made here, such as the complete absence of decoher- 
ence and symmetry breaking effects, i.e. essentially zero applied voltage, zero magnetic field, 
and zero temperature. The results are contrasted with those pertaining to the same system 
with the symmetries completely broken, but the specific form of the temperature or magnetic 
field dependence cannot be evaluated so simply (as in Ref. except for the case of a 

magnetic field in a multiply-connected geometry, with only simple Aharonov-Bohm phases. 
Also, only the leading-order, O(N), contribution to the conductance is evaluated here, in 
contrast with the transmission-matrix approach which often focuses attention on the 0(N°) 
corrections to the conductance and its fluctuations Another drastic simplification is 



the description of possibly complicated experimental geometries by a single scattering cavity, 
with the external leads considered as ideal. As already mentioned, work on extending the 
present approach to deal with more realistic situations is in progress [1^ . 



II. PLANAR DIAGRAMS FOR LARGE N UNITARY RANDOM MATRICES 

In this section we develop a diagrammatic technique which may be used to evaluate the 
average conductance of Eq. (0), and is accurate to leading order in 1/N. The diagrammatic 



S 



description of the unaveraged conductance is given in Fig. 2. Propagation of an electron 
through an ideal lead is denoted by a thin line with an arrow to the right; propagation of a 
hole is denoted by a line with a left-pointing arrow. An Andreev scattering event is denoted 
by a full semicircle. Each element of S carries two indices representing incoming and outgoing 
transverse modes in the leads, and corresponds to a dangling double line. These double 
lines will eventually be connected to each other in various ways, due to averaging over the 
ensemble of S. The fact that one of the indices carried by two consecutive S"-matrix elements 
is identical — the mode number is not changed by Andreev reflection — corresponds to the 
continuous single line which connects these two elements. The conductance itself, according 
to Eq. (|) and Eq. (j3|), corresponds to a sum over all possible "bubbles" with two electron 
lines leaving the left-most point and two hole lines arriving at the right-most point (see 
the figure), with an odd number of Andreev reflections occurring along the lower line and 
along the upper line. We use the convention that the upper line corresponds to the second 
T eh factor in Eq. (Q), and all expressions pertaining to it are thus complex-conjugated (or 
Hermitean-conj ugated) . 

In the first subsection below, we give a general discussion of averages of products of S- 
matrix elements, including our use of time-reversal symmetry to simplify the expressions, 
and the failure of a naive semiclassical approach to the evaluation of these averages. In 
the second subsection, we demonstrate that the corresponding averages may be described 
diagrammatically by introducing elements which couple k copies of S with k copies of S\ 
with weights as given in Fig. 3. Although the correlations represented by the higher-order 
couplings of Fig. 3 are obviously weaker than those represented by the k = 1, "simple" 
coupling, they affect more matrix elements, and thus contribute to the conductance of Eq. (0) 
in the leading order. It is remarkable that even though large random unitary matrices possess 
such more subtle correlations (when compared to the Gaussian ensembles of Hermitean 
matrices), an analogue of Wick's theorem still holds: a diagram composed of several of the 
elements of Fig. 3 carries a weight which, to leading order, is the product of the weights of 
the elements. 
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There is a close analogy to be drawn between the present diagrams and 't Hooft's analysis 
of quantum chromodynamics in the large N limit |26| : the double lines are analogous to gluon 
propagators and the single lines to those of quarks. This analogy is useful because the index 
structure is similar — each continuous single line carries an index which is to be summed 
over O(N) different modes (and each double line carries the corresponding pair of indices). 
Indeed, the basic rule that only "planar" diagrams contribute to the leading-order results 
carries over to the present situation: a coupling with a weight of order 1/N 2k ^ 1 contributes to 
the leading order only if it divides the plane of the diagram into 2k disconnected sub-regions, 
thus increasing the number of factors of N due to independent summations by 2k — 1. Note 
that the analogy with quantum chromodynamics is certainly incomplete — for example, the 
directions of the arrows in the diagrams is different, and there is no analogue here of gluon 
self-interactions (the couplings here connect S"-matrix elements, and have no dynamics of 
their own). Another important difference is the fact that the mode index here can belong 
to different leads, with Andreev reflections occurring only in the superconducting leads, and 
the external points (the left-most and right-most points) imposing restrictions to modes in 
normal leads — each index is thus summed only over a subset of the N possible values. 

As often done in diagrammatic quantum field theory, we will do the averaging in two 
steps. First, we restrict ourselves to connecting the double lines (S"-matrix elements) within 
T e h to each other (and similarly the double lines within T\ h to each other), i.e., we will sum 
over diagrams in which the upper and the lower lines in the conductance bubble do not 
interact with each other. Later, we will connect the double lines within T e h to those within 
Tg h , and form diagrams analogous to the diffuson contribution in the calculation of the 
conductance of disordered metals fl3fl . The generic features of these two steps, involving the 
"self-energy" and the "effective coupling", will be discussed in the third subsection below. 
Specific details which differ from one system to the next are deferred to the following section. 
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A. Generalities 



One of the simplest examples of random-matrix ensemble averages to be considered is 
that of a product of only two S-matrix elements: 

(<V4) = ^7 • ( 4 ) 

\ 3 3% l COE N + 1 

This is an exact result — valid for any N (the subscript COE implies that the system is 
time-reversal invariant). Summation over repeated indices is not implied in the present 
notation. In the following, it will be important to distinguish between two "leading-order" 
approximations of this result, the "naive" one being 



N+l N ' 1 ; 

and the other being 

1+^-1 () 

iV + 1 N v ; 

(note the use of different relational operators). Because each of the indices % and j is 
eventually to be summed independently over a set of a size of order N, the Sij term will 
turn out to contribute less by one factor of N, and is omitted in the "true" leading order. 
In the present work, the notion of "leading order" is always used in this second sense, unless 
the word "naive" is explicitly added. In this example, a term which is naively of leading 
order disappears in the eventual evaluation of the true leading order of the conductance. 
In contrast, below we will encounter terms which are omitted in the naive leading order, 
but eventually contribute to the leading order results because they contain less ^-function 
factors. 

For expressions involving arbitrarily many S-matrix elements, the naive leading-order 
results are easily found; for example they can be obtained from a semiclassical approximation 



[ 21 1 . In this approximation S is written as 
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where /i is an index enumerating classical orbits in the set {ij}, i.e. those that have initial 
conditions corresponding to mode j and final conditions corresponding to mode i (classical 
orbits which enter and leave the cavity at angles specified by j and i respectively, and 
propagate at the Fermi energy). The amplitude for propagation along the orbit \i is denoted 
by Ap, and its classical action by s^. 

Products of elements of S correspond to multiple sums over orbits, and averaging corre- 
sponds to dropping all terms which have non-trivial phase factors p7| . The approximation 
here is simply that % is small relative to the possible fluctuations in the actions s M , regardless 
of whether these fluctuations occur within a certain set {ij} due to the classically chaotic 
nature of the dynamics in the cavity, or between different realizations of a disordered poten- 
tial in the cavity. Consider, for example, the case without time-reversal symmetry, denoted 
by the subscript CUE (Circular Unitary Ensemble). In this case we have 

famous " ( £ A.Ale^A * ( £ \AA = i . (8) 

The last equality follows from the correspondence between |A M | 2 and classical probabilities, 
and the assumption of ergodicity, i.e. that the particle "forgets" its initial conditions once it 
enters the cavity, and is equally likely to leave the cavity in any outgoing mode or direction. 

For products of many elements of S and S\ the obvious generalization is to pair each 
element of S with an element of in all possible manners, so that cancelations of the 
actions can occur. For example, 

{s i3 s) k s kl st) cuE * + • (9) 

Such expressions may be obtained by joining the double lines in the diagrams using only 
the first element of Fig. 3, as displayed in the first few diagrams in Fig. 4 [specifically, that 
of Fig. 4(a), the first two in Fig. 4(b), and the first five in Fig. 4(c); in the last case the 
sixth possible permutation is omitted because it contributes at a higher order in 1/N]. This 
corresponds to ignoring the unitarity of the matrices, and taking the results for Hermitian 
matrices instead (the Gaussian Unitary Ensemble). As we will see in the next subsection, 
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an additional term of order 1/N 3 must be added to the right hand side of this equation in 
order not to violate unitarity, and such "naively" higher-order terms eventually affect the 
leading order of the conductance. 

In the case with time-reversal symmetry (COE), the action differences — s u vanish 
identically not only for v equal to \i G {ij}, but also for v equal to /i T G {ji}, which is the 
orbit related to fi by time reversal. This leads to an extra factor of 1 + 8ij in the second-to- 
last and last expressions in Eq. (§|), and similarly to additional terms in Eq. (|9|). However, 
the "true" leading order result is independent of the symmetry, e.g., 



The same holds also for averages involving arbitrarily many matrix elements alternating 
from S and S^. For the naive leading order contribution to the "true" leading order, this 
follows directly from the semiclassical approximation discussed above, as can be seen by 
imagining all possible ways of connecting all the matrix elements with simple double lines 
- whenever we use the symmetry of S to generate additional terms, we break the planarity 
of the diagram (the corresponding double line is crossed) and thus generate a contribution 
of a "true" higher order (containing additional 5-function factors). In the next subsection 
we will see that the terms which are of both naive and "true" leading order, together with 
the property of unitarity, suffice to determine the terms which are naively of higher order 
but contribute to the leading order of the conductance. Thus, as both the COE and the 
CUE are ensembles of unitary matrices, Eq. ( |TUD and its many-S'-and-S^ generalization 
must hold also beyond the naive leading order. Note that it is of course crucial to use 
and not S* in Eq. ([[(]). For example = 5ij/N and would contribute in a higher 



Below we will use the time-reversal symmetry property, S T = S and S* = S\ to rewrite 
all our expressions in terms of 5* and only (they naturally appear in alternating order), and 




(10) 
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between the amplitude for a hole produced by a single Andreev reflection and that produced 
by three consecutive Andreev reflections (with the multiplicative factors of i omitted). This 
term is first transformed to read (^S^SjiSjiSimS^SnkJ , and then the leading order in N 
result for its average is obtained within the CUE. Such averages over the unitary ensemble, 



i.e. over the unitary group U(N), have been performed, e.g., in Ref. PBfl , for general N. 
We proceed to follow a somewhat heuristic diagrammatic approach for evaluating them to 
leading order in 1/N, referring to this work for a mathematically rigorous derivation. 

If time-reversal symmetry is broken, the replacement of S* by is not permissible, 
and this affects the leading order results for the conductance, as already emphasized in 
the introduction. The symmetry may be broken by giving the electrons and holes a finite 
excitation energy, due to either a finite bias voltage or a finite temperature. In this case, 
the electrons travel with an energy Ep + e and the corresponding holes with an energy of 
Ep — e (where Ep denotes the Fermi energy and e the positive excitation energy). Thus S* 
is in fact the complex conjugate of S evaluated at a different energy, and the conductance 
may be obtained by assuming that S and S* are two independent members of the CUE, if 
e is greater than the range in energy over which correlations in S die out. This range or 
correlation energy, also called the Thouless energy, is easily evaluated from Eq. (0): one 
notes that the derivative of the classical action with respect to the energy is equal to t M , 
the duration of the orbit \x. All of these orbits are of lengths of the order of t esc , the escape 
time from the cavity, and thus the Thouless energy is equal to h/t esc . 

Another effect which breaks the symmetry between electrons and holes involves applying 



a magnetic field, stronger than the corresponding correlation field p9|| . Here S* is the 
complex conjugate of S (evaluated at the same energy), but is distinct from . However, 
the leading-order results for the conductance may still be obtained by considering them as 
independent members of the CUE (see, e.g., Ref. fT[|] ). This can easily be demonstrated using 
the diagrammatic language — all the planar diagrams have in this case the property that all 
the double lines connect S"s to S^'s and if>*'s to S T, s, with no intermixing [the argument is 
essentially the same as that used in the discussion of Eq. (|i0|) 1. In the physical discussion of 
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the following section, the results for the symmetric case are compared with those pertaining 
to such time-reversal non-symmetric situations, which turn out to be much easier to evaluate 
(cf. the last paragraph of this section). 



B. Generalized couplings for unitary matrices 

We thus would like to calculate expressions of the type ^Si a SljSjbSl k . . . Si c Sl^j, where 
(. . .) denotes the average over the unitary group U(N). We have used letters from different 
sections of the alphabet to emphasize the fact that, for elements of this group, the first 
(or upper) index and the second (or lower) index of S transform differently (covariantly 
and contravariantly) under U(N) (the reverse holds for S^). It follows directly from the 
invariance of the probability distribution of S under such U(N) transformations that the 
averaged product of any n elements of S and m elements of may be written as 

/ a a q at ct at \ _ 

^101^202 • • • °i n a n b 1 ,j 1 °b 2 j2 ■ ■ ■ °b m j m / ~ 

$nm Cu,v ^1^(1)^2^(2) • • • ^inju(n) ^1^(1)^12^(2) ' ' ' ^0„6 V („) J (11) 

u,v 

where u and v are permutations of n elements and C u>v are coefficients to be determined 
(the average vanishes of course for m different from n). Thus, ^-functions such as 5i a or 
products of them will never appear. 

These products of 5-functions may be represented diagrammatically using the elements 
of Fig. 3, with a simple k = 1 coupling for every 1-cycle in the combined permutation 
uv ~ 1 , a four-sided (k = 2) coupling for every 2-cycle, etc. The condition of planarity of 
the diagrams does not enter at this stage because we have not yet distinguished between 
contributions of different orders. The purpose of this subsection is to demonstrate that the 
leading-order expressions for the C u , u coefficients can be expressed as products of weights 
corresponding to the individual cycles in and that these weights are indeed those given 
in Fig. 3. Interestingly, we show that all the necessary properties of the C U; „'s follow from 
the unitarity of S, i.e. from the equalities J2 a SiaSlj = Sij and J2i SliSib = $ab- The unitarity 
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of S corresponds physically to the conservation of the number of particles in the cavity, i.e. 
to the fact that the incoming flux of electrons (or of holes) is equal to the outgoing flux. 

We will now calculate the required averages iteratively, starting from averages of only 
two elements and gradually building up. Start with 

(S m Sl) = 1 S |. (12) 

The equality here follows from the unitarity of S as we can see by summing both sides of 
the equality over a (or i). This is represented in figure 4(a) (note that in this figure the 
arrows have been omitted, and the matrix elements appear in counter-clockwise, rather than 
clockwise, order). We have introduced a quantity 72 defined to be 1 for later convenience. 
We next write 

/C Ct q 0t\ ~ jjj + 5ab I ^ 4 (lO\ 

\^ia^aj^jb^bi/ ~ N 2 N 3 

where the first term corresponds to the naive leading order result, and the second term is a 
possible correction of naively higher order. It involves 74, which is an unknown "coupling 
constant" of a four-sided interaction [see Fig. 4(b)]. To determine 74 we sum this equation 
over b. The left hand side gives 5ij (S ia Slj^ = ji, the right hand side ( 1+Si ^ + ^ , and thus 
we obtain 74 = — 1 = —7272- 

To see the recursion pattern, we go one step further and consider 

(SiaSljSjbS^SkcSl^ = ( ) + (5 ac + 5 ik + 5 bc + 5 jk + 5 ab + <%) + , (14) 

see Fig. 4(c). Here ( ) represents the naive leading order term — a sum of terms each 

involving two Kronecker deltas, for example, 5 a bSbc/N 3 or SikS a b/N 3 — which are generated 
by using only the "simple" 72 couplings. The second term, involving a single 5 function, is 
due to diagrams with one "simple" coupling, 72, and one four-sided vertex, 74. The third 
term corresponds to the six-sided vertex which must be introduced as well. 

One can easily see that the terms required by unitarity at (say) the c vertex are correctly 
reproduced by all the diagrams with a simple double line connecting Sk c with S^, and 
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all the previously derived elements (in this case 72 and 74) connecting the other four S- 
matrix elements. Of the twelve terms in Fig. 4(c), only three are of this type, and the 
other nine terms must thus cancel each other when summed over c. There are two types of 
cancelations between these "unwanted" terms. The first type occurs, for example, between 
the term 8 ab 5 bc /N 3 and the term — <5 ab /iV 4 . More generally, such cancelations occur due to 
the weights of the lower-order couplings already determined (in this case the weight 74), with 
a "distant" part of the diagram (in this case the simple coupling leading to the 5 ab factor) 
playing only a passive role. There are three pairs of diagrams in Fig. 4(c) which cancel in this 
way. The second type of cancelation accounts for the remaining three diagrams, including 
the one with the new six-sided vertex, whose weight 7 6 we need to determine. The diagram 
with the six-sided vertex involves no 5-functions, and the only terms which upon summing 
over c would produce no Kronecker delta are evidently the two terms 5 bc + 8 ca . Thus, we 
have the recursion relation 7 6 = —27274. 

Speaking picturesquely, we see that we have a combinatorial problem in which in general 
2k diners are seated around a round table and engaging in conversation. Note that the indices 

i,a,j, are to be associated with the spaces between the diners. Each diner is identified 

by two indices: thus, diner [ia], whom we may think of as a gentleman, is associated with 
Si a , diner [aj], whom we may think of as a lady, is associated with S^-, and so forth. The 
number of diners 2k is required to be always even. In the case 2k = 6 we just discussed, 7 6 
corresponds to all six diners engaging in one conversational group, 7472 corresponds to four 
of the diners engaging in one conversational group with the other two speaking only to each 
other. The rules of large N tells us that etiquette dictates that the "lines of conversation" 
cannot cross. 

Proceeding in this way, we see that for the general case the second type of cancelation 
will occur if 

fc-i 

12k + ^2l2ll2{k-l) = . (15) 
1=1 

We can readily interpret the terms as follows: in the first term all of the 2k diners are 
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engaged in one conversation, while in the second term the group has broken up in two, with 
the two diners who share the index on which we are summing (which was index c in the above 
example) belonging to different groups. The sum is over all terms with a single 5 function 
connecting that index (c) with all others of its type (covariant or contravariant). There are 
always an even number of diners in any conversational group, the number of "males" being 
equal to the number of "females" — the subscript of 7 is always even. 

This last fact clearly suggests defining Cj = (—l) J+r y2j, in terms of which we have the 
recursion relation 

k-i 

Ck = C 3 C k-j J ( 16 ) 

J'=l 

with the explicit solution 

_ (2k)\ 

[another way of stating this result is: = (4 — 6/k)ck-i, C\ = 1]. Actually, for our later 
applications we don't need the explicit solution for as much as the generating function 
defined by w(x) = Y^T=i c k x ■ The recursion relation above immediately implies 

w = x + w 2 ; (18) 

the extra term x on the right hand side can be seen by noting that the series for w(x) starts 
with x for small x and hence w 2 starts with x 2 . We thus directly obtain 



w 



(l - VT^) /2 , (19) 



which implies fllTD - In the appendix, we give an alternative derivation, which involves gener- 
alizing the energy-dependent Green's functions often used for random Hermitean matrices 
to the case of unitary matrices. 

In the above, we have tacitly assumed that the weights of diagrams with several coupling 
elements are given by the product of the weights of the elements (those of Fig. 3). Our 
argument so far does not demonstrate that this is the only way to maintain unitarity. 



However, the unitarity relations between diagrams with 2k matrix elements and with 2k — 2 
matrix elements do in fact suffice to determine the weights uniquely ||28|| , and therefore the 
weights we have found are the correct weights. To see this, imagine that all weights of all 
diagrams with 2k — 2 or less matrix elements have been determined, and that you need to 
determine the weight of the diagrams with 2k elements. Start from all diagrams having 
at least one simple double line connecting two adjacent "diners" or matrix elements. By 
summing over the index shared by these two diners, one immediately gets the weight of that 
diagram as 1/N times the weight of the corresponding diagram with that double line omitted 
(as well as the two diners it connects). For diagrams which do not have such pairs, one can 
start from a foursome of adjacent diners, and sum over any one of the indices internal to 
that foursome. Obviously this can be continued until all diagrams of 2k elements have been 
accounted for, with the last diagram serving to determine according to Eq. (|T5|). Note 
that the actual C u>v for general N have more complicated expressions, with special cases for 
n > N. However, in the large N limit, corrections to the above analysis (e.g. due to non- 
planar diagrams) are smaller than the leading term by a factor of 1/N 2 . This high degree 
of accuracy of the leading order is specific to the case without time-reversal symmetry - 
corrections of relative order 1/N would be present in the COE. 



C. The self energy S and the effective coupling T 

According to the diagrammatic rules established above, only planar diagrams contribute 
to the leading order results. In the present subsection we use these rules to define a "self- 
energy" and an "effective coupling" in a manner analogous to similar quantities which appear 
in diagrammatic analyses of perturbative field theory. We thus group together diagrams 
which differ from each other in their internal structure, but not in the way they connect 
to the rest of the diagram (or the rest of the "plane"). There are two relevant families 
- the "one-particle irreducible" diagrams are summed in Fig. 5(a), and the "two-particle 
irreducible" diagrams are summed in Fig. 5(d). These elements may then be used to calculate 
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the "averaged Green's function" G, the "diffuson" D, and eventually the conductance a (see 
Fig. 5). Note that although the analogy to the usual perturbation theory of disordered 
metals is both striking and useful, there are also important differences: the scattering is 
included in a very different way, and unlike a typical Green's function in quantum field 
theory, there is no energy variable involved in G or E (at least at the present level). 

Following from Fig. 5, we have a "Schwinger-Dyson" equation for the Green's function 

G = G + G SG , (20) 

or equivalent ly 

G- 1 = G 1 - E (21) 

(from here on we use = rather than =, although it should be remembered that results 
are being evaluated only to leading order). The notation here is such that Gq represents 



propagation with a single Andreev reflection p0[ — it is a two-by-two, off-diagonal matrix 



in electron-hole space, whereas it is diagonal in mode-space. Similarly, S and G are off- 
diagonal in electron-hole space. Note that the self-energy £ differs from the other quantities 
in this equation in that it is proportional to the unit matrix in mode space — the diagrams 
of Fig. 5(a) are independent of the (single) mode index carried by the external lines. 

It is convenient to parameterize £ by a single complex parameter /, 

/0 if*\ 

£ = , (22) 

w o ; 

where the structure in mode space is not explicitly written. The fact that the electron-hole 
element if is minus the complex conjugate of the hole-electron element if* follows directly 
from the fact that a single Andreev reflection, and thus Go and G, also share this structure. 
One may refer to / as the amplitude for Andreev reflection in the cavity, because following 
from Fig. 5(a), to leading order in N and for the averaged behavior, an electron entering the 
cavity from any one of the leads may be considered to be Andreev reflected directly into the 
time-reversed hole trajectory, with a probability amplitude / (and the extra phase factor 
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The Green's functions Go and G have a more complicated internal structure. To maintain 
the generality of the present section, we consider here a cavity connected to several supercon- 
ducting and normal leads. The bare Green's function, Go, is parameterized as in Eq. (p2|). 
with / replaced by e - ^, using different values of (f) for modes in different superconducting 
leads, and / replaced by (and hence Go = 0) for modes in normal leads. The corresponding 
elements of G are, according to Eq. (pT|) , similarly described by a parameter g = l/(e 1 ^ + f*), 
and since <fi varies from one superconducting lead to the next, so does g. This total ampli- 
tude for transforming an electron into a hole can be understood as a sum over terms with 
an odd number of Andreev reflections: ig = ie~ 1 ^ + i^e^ j*e~ 1 ^ + i^ e - l <t>f* e -i<t>j* e -^ _|_ 

The specific values of g, and the number of modes for which each one of them occurs, 
depend on the geometrical structure of the system. Such specifics will be discussed in the 
following section. However, one may make progress by introducing a notation for the trace 
of each of the off-diagonal blocks of G, which reads iNa = J2j Gj,h;j,e, an d leads directly to 

a = E^n = E ^r^^p (23) 

where the sum over n runs over the superconducting lead(s). The trace of the hole-electron 
block is similarly equal to iNa*. We proceed to derive a general relationship between a 
and /, which taken together with Eq. (|23|) will eventually allow us to find / for particular 
systems. 

According to Fig. 5(a), the "self energy" E is given by a sum of terms involving higher 
and higher order couplings. Each r y 2 k coupling divides the plane of the diagram into an 
external part and 2k — 1 internal "triangles" , which give rise to factors of a and a* (traces 
over Ghe and G e h). Mathematically, this relationship reads 

/ = 7 2 a - 7 4 a|a| 2 + 7ea|a| 4 + • • • = E 72fc" (-|a| 2 J , (24) 

fc=i 

or, in terms of the c^s, 

/ = ]Tc fc aM 2fc - 2 . (25) 
fc=i 
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Comparing with (|T^) we find 



/ = ^(i-v^«) 



(26) 



Actually, a more useful form comes directly from Eq. ( |T8| ) with the substitution w — > a* f 
and x — > \a\ 2 (using the fact that the complex phase of / is equal to that of a): 



a result which we will use repeatedly. 

We now turn to a similar consideration of the "two-particle irreducible" diagrams, i.e. 
the "effective coupling" T depicted in Fig. 5(d), and its relationship with the "diffuson" 
D of Fig. 5(e). The structure in mode space is again trivial, with both V and D having 
two independent mode indices on the left and right (in fact one should consider V and D 
to have four mode indices, with the right two and the left two being equal to each other). 
However, the structure in electron-hole space is more complicated than before, having both 
diagonal and off-diagonal elements. Before any summation over the external mode indices, 
the diagrams in Fig. 5(d) and Fig. 5(e) are of order 1/N, and we define T and D to be N 
times the corresponding diagrams. We parameterize T by T = (T^I + T c t), where I (often 



Here (for "diagonal" ) is the probability for an electron to remain an electron, and a hole 
to remain a hole, while T c (for "change") is the probability for an electron to become a 
hole and vice versa. Similarly, we write D = (Ddl + D c t) . According to the diagrams of 
Fig. 5(e), we have 



omitted) 




denote the identity matrix and the first Pauli matrix respectively. 



D 



Y + YtD 



(28) 



or 



D 



i 




(29) 



where 
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' = E^k|V (30) 

n 

represents the two Green's functions G and G* connecting T to D. The J2 n W n appears 
explicitly in our two-by-two matrix notation, whereas it would be implicit in the matrix 
multiplication in Eq. (^) if multiplication in the mode sub-space were implied [the factor 
of 1/N is introduced in t to account for our definition of T and D as 0(N°) objects]. Again, 
the specific values for t and D depend very much on the geometric structure, and are left 
for the following section. Here we proceed to derive the general relationship between V and 
a (and therefore /), which, as in the case of the self-energy S, is the only aspect of G which 
enters in the definition of V. 

Let us then look at the graphs contributing to T, Fig. 5(d), in particular the terms with 
2k double lines joined by a coupling constant j2k- Those graphs in which an odd number of 
double lines land on the top line (and thus necessarily with an odd number of double lines 
landing on the bottom line) contribute to T^, and not to r c . There are k such possibilities 
(one double line can land on the top line, three double lines, and so on, up to 2k — 1 double 
lines). In contrast, those graphs in which an even number of double lines land on the top 
and bottom lines contribute to r c . There are k — 1 such possibilities. We thus find that 

oo 

T d = Y,Ck\a\ 2k ~ 2 k, (31) 

k=l 

where in each contribution to the signs in (— \a\ 2 ^j cancel with the sign of 72^- Com- 



paring with (^) we see that T d = ,(} , and thus 



I 1 _ i + l/l 2 (32) 

'1-4 a 2 1 — 1/1 



from Eqs. (|26l) and ( p7|) . Similarly, we have 

00 

-T c = Y,c k \a\ 2k ~ 2 (k-l). (33) 

k=l 

Note the minus sign in the last equation. To see its origin, consider the term involving 
four double lines with two landing on the top single line and two landing on the bot- 
tom line: it is given by 74(20) (—ia*) for the electron-hole term, and 74(20;*) (—ia) for the 
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hole-electron term (they are equal; recall that expressions pertaining to the top line are 
complex-conjugated). Due to the negative sign of 74, this gives a negative contribution to 
T c . Next, consider a typical term involving six double lines, say with two double lines going 
upwards and four double lines going downwards. The contribution to T c is given by, say, 
7 6 (i 3 a;|a;| 2 )(— ia*), but 7 6 is positive. Proceeding in this way, we see that T c is given by 
a series with negative coefficients. In contrast, Td is given by a series with positive coeffi- 
cients. Consider for example the term just mentioned with six double lines. To obtain a 
contribution to we have to move, for example, one double line from the bottom to the 
top part of the diagram. In the expression just given, one factor of ia is thus changed to 
— ia. As already emphasized in the introduction, various signs of this type, and the factors 
of i which generate them, are of crucial importance to obtaining our final results, and reflect 
fundamental aspects of the physics involved. 

To continue with our evaluation of r c we obtain by its defining series (|33| ) that 



r r / m 2 (i + m 2 ) m 

Putting Trf and T c together and evaluating T~ l , we find the remarkably simple formula 

^ = ^ = rn7FF (1+l/W - (35) 



As we saw from (|29|) , we need the inverse of T, rather than T itself. 



In the following section, we simply apply these formulas — Eqs. fl23|) , (|2~7[), fl29|), (|30| 



and (|35D — to a specific geometry, and proceed to evaluate diagrammatically the average 
conductance of Eq. (||). 

Before continuing, we mention how the procedure is changed when one wishes to evaluate 
the conductance in the absence of time-reversal symmetry. According to the discussion at 
the end of subsection A above, in this case the leading order contributions cannot contain 
couplings connecting electron scattering elements and hole scattering elements along the 
bottom line (or the top line) of the diagrams. Thus there are no contributions to X (to 
leading order), and we find that / = in this case (and thus also a = 0). Similarly, only the 
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first "simple" diagram in Fig. 5(d) is a legitimate contribution, and T = = /. Lastly, 
we have G = G in this case, and \g n \ = \e~ t</>n \ = 1 should be used in Eq. (p0|). In order to 
obtain the conductance, one still needs to invert Eq. (p5| ) as in the symmetric case. 



III. A SIMPLE APPLICATION 

The simplest possible dirty mesoscopic N-S system arguably consists of a disordered 
cavity or junction with just one normal and one superconducting lead. In such a system 
the phase of the superconducting order parameter can be gauged away, and so we assume 
that = and / is real. We denote the relative width of the superconducting lead by 
n s = W\/N, so that the conductance we seek is a simple function of < n s < 1, multiplied 
by the appropriate units and the total number of modes N. 

From the Schwinger-Dyson equation, Eq. (^1|), we see that 

9 = T~V~e > ( 36 ) 



or equivalently 



1 + / 

from Eq. (p3|). Putting Eqs. ( pT|) and (|36|) together, we obtain 



1 + / 



« = 7^7 , (37) 



1 + P 1 + / ' 
and thus we determine / completely in terms of n s : 



(38) 



1 + 4n s - Anl - 1 



/ = ; ^ • (39) 

2(1 -n a ) 1 ; 

The plus root is chosen so that as the number of modes n s in the lead connecting the 
cavity and the superconductor goes to zero, / should vanish as is physically reasonable. The 
amplitude / grows with n s , and approaches 1 when n s approaches 1. In particular, for the 
symmetric case, n s = | we have / = \[2 — 1. 
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We proceed to evaluate the diffuson for this system. Using Eq. fl29|) with 



* = n£ ' T = (i+/4+/) r ■ (40) 

we find through simple arithmetic that 

(1+ /»)(! + /)" / f(l-f) \ ,,,, 

D= d + 2/-/') l / + 7TT7r T j ' (41) 

Now we are ready to compute the conductance, as represented diagrammatically in 
Fig. 5(f). As shown in the figure, there are two distinct terms which contribute to this 
quantity — one in which the electrons are Andreev reflected "individually" (which involves 
only E), and one in which the top line and the bottom line in the diagram are connected, 
which involves the diffuson D. 

In the case of the non-diffuson contribution, one must note that the appropriate terms of 
G must contain only a single copy of S, as opposed to the infinite series reflected in Eq. (|36|). 
This is due to the fact that the external indices represent modes in the normal lead, and 
cannot be directly Andreev reflected (Go for the normal lead vanishes — excitations leaving 
the cavity through it never return). Thus, the non-diffuson contribution to the conductance 
is given by 

2 2 

a E = A e -N{\ - n s )f = 4jN(n s - /) . (42) 

This is plotted in Fig. 6 function of n*. 

The conductance due to the diffuson is given by 

°d = *jN(l - n s fD c = AjN{1 - n s ) 2 /I 1 ' P f2 , (43) 

which is also plotted in Fig. 6. Note that to obtain the non-diffuson contribution to the 
conductance we multiply by the number of modes on the external lead W<i = (1 — n s )N, 
while to obtain the diffuson contribution we multiply by the square of this number. In the 
non-diffuson contribution to the conductance, an electron injected into a given mode comes 
back as a hole in the same mode: this is known as the "giant backscattering peak" |IU |. In 



contrast, in the diffuson contribution, the hole comes back in general in some other mode. 
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The total conductance is of course given by the sum of Eqs. (f42|) and (f43l). We find that 
when we add these two contributions to the conductance, we obtain the surprisingly simple 
expression 

ctns = W + (Td = (1 - 1 ] . (44) 

h \ Jl + 4n s - 4n? } 



Indeed, if we introduce the asymmetry parameter a by n s = (1 + a)/2, so that a = for a 
symmetric cavity, we obtain simply 

a NS = (l - J—z) ■ (45) 



a 



This is again plotted in Fig. 6. For the case of equal widths of the leads, a = 0, this has 



been derived using the transmission-eigenvalue approach in Ref. |31| . Incidentally, note 
that <jns is symmetric with respect to interchanging the widths of the leads, a <-> —a. In 
the transmission-eigenvalue approach, this follows directly from the fact that (7ns depends 
only on the transmission eigenvalues |32| (even when these are appropriately redefined for 
rectangular transmission matrices). 

For comparison, we consider the case in which the charge conjugation symmetry between 
electrons and holes is broken, by a finite bias voltage V, a finite temperature, or a finite 
magnetic field. As explained above (see the end of Sec. II A and Sec. II C) we may then 
denote the scattering matrix of holes by a unitary matrix 5", which to leading order can 
be taken as unrelated to S, and average over S and S' separately. The situation becomes 
enormously simpler: all the non-diffuson graphs cease to exist, S = 0, and we need only 
to evaluate the diffuson contribution with T = I and G = Gq. We are left with the ladder 
graphs shown in Fig. 7, with the number of rungs on the ladder restricted to be even. The 
conductance is thus given by (1 — n s ) 2 (n s + + n 5 s + ...) = n s (j^ L ) = 2^HiQ ' wn i c h we 
plot in Fig. 8. Notice that it is not symmetric in the asymmetry parameter a. We also plot 
the difference between the total conductance with and without charge conjugation (or time 
reversal) symmetry. We see that for n s < |(1 + a c ) ~ 0.65 charge conjugation symmetry 
actually lowers the conductance (a c = 0.2955 ... is the real solution of a? c + a 2 c + 3a c — 1 = 0). 
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Incidentally, the total conductance without charge conjugation symmetry can be calcu- 
lated easily using elementary physics. Since the scattering of the hole in the cavity is no 
longer strongly correlated with the scattering of the electron, the system is analogous to 
two independent cavities connected in series, the motion of electrons being represented as 
excitations in the first cavity, and holes being represented as excitations in the second. The 
conductance of such a two-cavity system is given, to leading order, simply by the series ad- 
dition of the resistances corresponding to the leads that connect them to each other and to 
the "electron electrode" and its copy, the "hole electrode": a 2C = 2^(J¥ 2 ~ 1 + WT 1 + W2 1 )' 1 
(the subscript 2C implies 2 cavities). According to Eq. @, the conductance of the N-S 
structure is given by twice this value, which is, as noted above, 

non— symmetric a at s t Ar\ 

= 4iV-n s — . (46) 

Another comparison may be made with the conductance the same structure would have 
if both electrodes were normal, = (2e 2 /h) Nn s (l — n s ). However, experimentally it is 
usually much easier to observe the crossover from a^s to a ^° -s y mmctnc than that from a^s 

tO <TjV- 

The results reported here are in complete agreement with previous theoretical work. For 
example, the fact that breaking electron-hole symmetry may either decrease or increase the 



conductance, ctjvs, is known, at least in principle [Q. This situation suggests the follow- 
ing experimental challenge: to fabricate a device with an ergodic scatterer connected to 
electrodes through leads with tunable widths, for example by changing a gate voltage in a 
two-dimensional electron gas device, such that ctns would exhibit positive magnetoconduc- 
tance for some gate voltages and negative magnetoconductance for other gate voltages. The 
challenge here is considerable because fabricating good contacts between a two-dimensional 
electron gas and a superconducting electrode is not easy — there are invariably large mis- 
matches in the Fermi velocity, not to mention Schottky barriers, etc. In fact, the goal of 
fabricating such contacts has been one of the driving forces behind the development of this 
field in the past few years. Note that related normal-normal conductances have already been 
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shown to exhibit changes in the sign of the magnetoconductance as a function of geometric 
parameters of the system, in Ref. [[J. 

It is interesting to compare our results with the experimental data of Kastalsky et al. 
]1|], who have studied an N-S system in which a relatively large piece of "normal metal" 
fin this degenerate semiconductor) was in contact with a superconductor and a 

normal reservoir (the semiconducting substrate). The conductance of this junction a^s was 
observed to grow by a factor of about 1.8 upon decreasing an applied magnetic field at 
low temperatures (the temperature dependence is complicated by not having a very large 
separation between the Thouless energy and the superconducting gap). This compares 
favorably with the results shown in Fig. 8 for large n s — in this case / approaches 1 while 
the resistance of the wide lead, Wi, becomes negligible; the conductance due to "direct" 
Andreev reflections (os) is then equal to twice that which is obtained when the holes must 
"find their way to the external electrodes" independently of the motion of the electrons 
that produced them [expanding Eqs. ( fP|) and ( pj|) near n s — 1 gives a^s — > and 
<t ns ~ symmctnc — > ^-W2 respectively]. This experiment was in fact one of the first available 
in the subfield of mesoscopic N-S structures discussed here, and was initially very hard 
to interpret. The interpretation which was eventually accepted 0] put forward precisely 
the ideas of multiple Andreev reflection which are formalized in the present work and its 
predecessors. 



IV. CONCLUSIONS AND OUTLOOK 

We have developed a new technique for evaluating the dissipative conductances of disor- 
dered mesoscopic systems attached to normal and to superconducting leads. It is based on 
a generalization of the planar-diagrammatic technique of large-A random matrix models, 
which allows one to consider unitary (as opposed to Hermitean) matrices, specifically the 
scattering matrix of a disordered mesoscopic grain. The leading order results for the con- 
ductance are affected by mesoscopic coherence, because of the possibility and importance of 
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an exact symmetry between electrons and holes at the Fermi level. 

In the present work, we have demonstrated the method only by application to the simplest 
possible disordered mesoscopic N-S system: a scattering cavity attached through ideal leads 
to one superconducting electrode and one normal electrode. This calculation can only serve 
as a "toy-model" description of actual experimental geometries; it is necessary to consider 
many possible complications, such as the effects of potential barriers, in order to check, e.g., 
whether there is any physics behind the agreement with the experiment of Ref. [|IJ which was 
mentioned at the end of the previous section. This underlines the importance of generalizing 
the present method to deal with more complicated geometries (in analogy with the "circuit 
theory of Andreev conductance" put forward by Nazarov |7|). Such a generalization turns 



out to be possible, and the results will be reported in future work |L6 . 

An additional possible direction for future pursuit is to evaluate the fluctuations and the 
next order corrections to the conductances. Specifically, the much-discussed 0(N°) term can 
in principle be found from similar large- N diagrams which break the condition of planarity 
exactly once. A third direction is to consider excitations with finite energies, or systems in 
a finite magnetic field. Such calculations could turn out to be much more difficult, because 
the scattering matrices depend on energy and magnetic field in a continuous manner, and 
it should be necessary to invoke ensembles of such two-parameter families of scattering 
matrices. However, the experience in mesoscopic phenomena may lead to a hope that the 



generic behavior in such ensembles would be simple and universal [35[], perhaps describable 



by simple modifications of the couplings used in the present method. 



APPENDIX 



Here we give an alternative derivation of fllTD - As customary in the literature on random 
matrix theory, we define the function of a complex variable z called the "one-point" Green's 
function by 

G ' 2 » = (h^T^hm) (47) 
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(which should not be confused with the quantity G in the main text). As before, (. . .) 
denotes averaging the N by N unitary matrix S over the CUE. Diagonalizing the unitary 
matrix S and denoting its eigenvalues by e %dj , j = 1,2, ....N, we may use the obvious fact 
that the eigenvalues are distributed uniformly over the unit circle to find that 

G(z) = f d6 1 = f J_M1 = 1 (48) 
1 ' J 2tt z - 2cos6 ^ z 2 ^ 1 (k\) 2 ' ' 1 ' 

where the sum converges for \z\ > 2. 

On the other hand, if we evaluate ( f4T|) by a diagrammatic expansion as shown in Fig. 9, 
we have the two equations 

G « = 7=W) < 49) 

and 

oo 

E(z) = 2j2l2kG 2k ~\z) (50) 
fe=i 

(the latter holds at least to leading order in 1/N). We recognize as just the Schwinger- 
Dyson equation again: here the bare Green's function Gq(z) is simply equal to 1/ z. In (50) 



we have used the definition of r )2k as the coupling constant involving 2k matrix elements. 
Note also the factor of 2: the first matrix element in the diagram for E, reading from left 
to right say, may represent either S or S^. This feature does not appear in discussions of 
Hermitean random matrices. 

Our goal here is to determine the coupling constants 7 2 fc, assuming that averages over 
elements of S and S< may indeed be written in terms of such diagrams, with a generalization 
of Wick's theorem implying a multiplicative property of the weights of diagrams involving 
more than one coupling (those of Fig. 3). We do this by eliminating z between fl4lf ) and 
(P3) and thus solving for E(z) in terms of G(z): 



Jl + 4G 2 (z)-l 
E(2) - ' G(z) ' < 51 > 

which is identical to Eq. ( fOP if we identify x = —G 2 (z) and w = —H(z)G(z)/2 [cf. also 

Eq. (P7|)]. Expanding the right-hand side as a series in G and comparing with fl5H| ) we 

obtain immediately j2k and hence in agreement with ([l] 
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FIGURES 



Fig. 1: Sketch of a simple normal-superconducting mesoscopic system, consisting of 
a single "grain" which is treated as an ideally ergodic scattering cavity, connected to a 
normal electrode (particle reservoir) and to a superconducting electrode through ideal leads 
of different widths. 

Fig. 2: An example of a diagram contributing to the unaveraged conductance of Eq. (§). 

Fig. 3: The different couplings, and their weights [see Eq. ( |i~7l) below]. These couplings 
are used to connect the double lines of Fig. 2. 

Fig. 4: Diagrams corresponding to (a) (Si a Sl^, (b) (SiaSljSjbSl^, and (c) 

/ n ct c ct c ct \ 
\ lD ia>-> a j>3jbO bk Dkc>->ci / ■ 

Fig. 5: Diagrammatic representation of: (a) the "self-energy" S, (b) the "bare Green's 
function" Gq, (c) the "averaged Green's function" G, (d) the "effective coupling constant" 
T, (e) the "diffuson" D, and (f) the conductance a of Eq. (Q). 

Fig. 6: The conductance of the simplest normal-superconducting system considered, in 
units of 4Ne 2 /h, as a function of the fraction of modes belonging to the superconducting 
lead (full line). The non-diffuson (dashed) and the diffuson (dotted) contributions are also 
shown separately, the former corresponding to "direct" Andreev reflections from the cavity, 
or the so-called giant backscattering peak. 

Fig. 7: The simplified diagrams contributing to the conductance when the symmetry 
between electrons and holes is broken. Only vertical simple double lines can be used, because 
the diagram must be planar and the matrices on the even rungs of the ladder cannot be 
connected to the matrices on the odd rungs. 
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Fig. 8: A comparison of the conductance of Fig. 6 (full line) to that obtained in the 
absence of electron-hole symmetry (dashed line). Also shown is the difference between the 
two, on an expanded scale (top panel). 

Fig. 9: Diagrams for the one-point Green's function, G(z) (thick lines), and its associated 
self-energy E(z) (shaded semicircles) of the appendix. The thin lines represent G (z) = 1/z, 
and carry no arrows. The double lines represent an element of S if the arrow is pointing 
away from the single line, and an element of in the opposite case. Although the quantities 
G(z) and £(z) are not to be confused with G and £ of the main text, the couplings are the 
same (see Fig. 3). 
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figure 1, Argaman and Zee 




figure 2, Argaman and Zee 




figure 3, Argaman and Zee 




figure 4, Argaman and Zee 




figure 6, Argaman & Zee 
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figure 7, Argaman and Zee 




figure 8, Argaman and Zee 



